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ABSTRACT 

We report on ii'-band interferometric observations of the bright, rapidly ro- 
tating star Regulus (type B7 V) made with the CHARA Array on Mount Wilson, 
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California. Through a combination of interferometric and spectroscopic measure- 
ments, we have determined for Regulus the equatorial and polar diameters and 
temperatures, the rotational velocity and period, the inclination and position 
angle of the spin axis, and the gravity darkening coefficient. These first results 
from the CHARA Array provide the first interferometric measurement of gravity 
darkening in a rapidly rotating star and represent the first detection of gravity 
darkening in a star that is not a member of an eclipsing binary system. 

Subject headings: stars: fundamental parameters — stars: individual (Alpha 
Leo, Regulus) — stars: rotation — infrared: stars — techniques: interferometric 



1. Introduction 

1.1. The CHARA Array 

Georgia State University's Center for High Angular Resolution Astronomy (CHARA) 
operates an optical/IR interferometric array on the grounds of Mount Wilson Observatory 
in the San Gabriel Mountains of southern California. The six light-collecting telescopes 
of the CHARA Array, each of 1-m aperture, are distributed in a Y-shaped configuration 
providing 15 baselines ranging from 34.1 to 330.7 m. Three of these baselines formed by 
the outer telescopes of the Y are in excess of 300 m. Ground breaking for the facility 
occurred in 1996 July, and the sixth and final telescope became fully integrated into the 
Array in 2003 December, thus signaling completion of construction and readiness for on- 
going science observations. Numerous technical reports and updates have been published 
throughout the design and construction phases of this project, the most recent of which is by 
McAlister et al. (2004), with many internal project reports available at CHARA's website^. 
A companion paper (ten Brummelaar et al. 2004) describes the full complement of technical 
and performance details of the instrument. 

With its long baselines, currently the longest operational K-hand baselines in the world, 
and its moderate apertures, the CHARA Array is well positioned for observations of main 
sequence stars not previously accessible to such high-resolution scrutiny. CHARA inaugu- 
rated its first diverse observing season in the spring of 2004 with a fully scheduled semester of 
science observations targeting stellar diameters, young stellar objects, and rapidly rotating 
stars. That semester was immediately preceded by an observing campaign directed entirely 
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at Regulus, a very rapidly rotating star expected to exhibit marked rotational oblateness and 
a resulting surface flux variation associated with gravity darkening. This paper is a report of 
the outcome of those observations of Regulus and represents the first scientific results from 
the CHARA Array. 

1.2. Regulus 

The star Regulus {a Leo, HR 3982, HD 87901) is of spectral type B7 V (Johnson & 
Morgan 1953) or B8 IVn (Gray ct al. 2003), and is a well-known rapid rotator. The difficulty 
of determining the projected rotational velocity Vsini is exemplified by the wide range of 
values published for this star, with velocities in the literature ranging from 249 ± 9 km s~^ 
(Stoeckley Carroll, & Miller 1984) to 350 ±25 km s"^ (Slcttcbak 1963). The first attempts 
to model realistically a line profile derived from a rotating star were by Elvcy (1930) who 
employed the method proposed by Shajn & Struve (1929) based upon dividing a photosphere 
into strips whose Doppler shifted contributions are weighted by the relative areas of these 
strips. Slettebak (1949) very significantly extended this analysis for rapidly rotating O and 
B stars to include the effects of limb darkening, gravity darkening based on Roche models 
incorporating the von Zeipel Effect (von Zeipel 1924), and differential rotation. Slettebak 
(1966) subsequently showed that the envelope of the most rapidly rotating stars of a given 
spectral type most closely approached the threshold of breakup velocities for stars of carly- 
to mid-B type. Regulus clearly falls in the domain of stars rotating close to breakup speed 
and can thus be expected to show significant rotationally induced oblateness and gravity 
darkening. 

The Washington Double Star Catalog^ lists three companions to Regulus, none of which 
has exhibited significant orbital motion since its first detection. The B component, ~ 175 arc- 
sec from Regulus, is HD 87884, a K2 V star whose proper motion, radial velocity and spec- 
troscopic parallax are consistent with its being a physical rather than optical companion. 
On the other hand a comparison of ages of components A and B has found the discordant 
values of 150 and 50 Myr, respectively (Gerbaldi, Faraggiana, & Balin 2001). Component 
B is accompanied by a faint companion {V = +13.1) comprising the subsystem BC that 
closed from 4.0 arcsec to 2.5 arcscc during the period 1867 to 1943 and has apparently not 
been measured in the last 60 years. At the distance of Regulus, the apparent magnitude of 
component C implies a star of approximate spectral type M4 V. Only a single measurement 
exists for a D component, more than 200 arcsec from Regulus, whose existence as a physical 
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member of the a Leo system might therefore be suspect. 

The Bright Star Catalog (HofHeit 1982) notes that Regulus is a spectroscopic binary. 
We have found no reference to this in the hterature except a flag in the bibliographic catalog 
of Abt & Biggs (1972) indicating the "progressive change" of radial velocity as reported by 
Maunder (1892). After conversion to modern units (Maunder used German geographical 
miles per second), the measured velocities declined from +40 km s~^ to —9 km s~^ between 
1875 and 1890. In light of the difficulty of accurately measuring velocities from hydrogen lines 
compounded by the very high rotational speed and the state of the art of the technique more 
than a century ago, the conclusion that those data indicate orbital motion seems unlikely. 
Furthermore, the General Catalogue of Stellar Radial Velocities (Wilson 1963) does not flag 
the star as a spectroscopic binary. On the other hand, Regulus may represent a case of 
neglect by radial velocity observers due to its brightness and extreme rotational broadening. 
If the velocities were to hint at a closer companion than the AB system, the corresponding 
orbital period would be on the order of a year with a separation in the range 50 to 70 mas, the 
regime of detectability by speckle intcrfcrometry. However, no closer companions have been 
reported from lunar occultation observations or from the CHARA speckle inter feromctry 
survey of bright stars (McAlister et al. 1993). In this angular separation regime, the fringe 
packets from Regulus and a very much fainter close companion would be non-overlapping 
in the fringe scan, and the companion's presence would have no effect on our analysis of 
stellar shape and surface brightness. One of the goals of this study has been to determine 
the orientation of the spin axis of Regulus. Owing to the lack of sufficient motion within the 
Regulus multiple star system and the apparent lack of any unknown closer companions, it 
will not be possible to compare rotational and orbital angular momentum vectors. 

Regulus possesses an intrinsic brightness and relative nearness to the Sun to make it an 
ideal candidate for diameter determination from long-baseline interferometry. Therefore, the 
star was among the classic sample of bright stars whose diameters were measured by Han- 
bury Brown and his colleagues at the Narrabri Intensity Interferometer nearly four decades 
ago (Hanbury Brown 1968; Hanbury Brown et al. 1967; Hanbury Brown, Davis, & Allen 
1974). In the first of these papers, the implications of a very high rotational velocity for the 
derived diameter were clearly recognized and discussed semi-quantitatively, but the result- 
ing diameter of 1.32±0.06 mas, measured at A = 438.5 nm for a uniform disk (UD), was 
judged to incorporate insufficient baseline sampling to be interpreted as other than circularly 
symmetric. A value of 1.37±0.06 mas was determined for a hmb-darkened disk (LD). Be- 
cause gravity darkening brightens the poles and thus counteracts the star's oblateness, it was 
thought that the maximum change in apparent angular diameter within the primary lobe 
was only 6%, which was judged to be marginal to detect in those observations. Furthermore, 
measurements beyond the first lobe were impossible because of the limited sensitivity of the 
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Narrabri instrument (in which the response varies as the square of the visibihty) . 

The other direct measurements of the angular diameter of Regulus come from lunar 
occultations. Although an early measurement (Berg 1970) yielded a diameter of 1.7±0.5 
mas, the only two reliable measurements of Regulus have been made by Radick (1981) 
and Ridgway et al. (1982). Radick (1981) obtained a diameter of 1.32±0.12 mas (UD) at 
A = 576.8 nm and 1.37±0.11 (UD) at A = 435.6 nm. These translate to 1.36 and 1.43 
mas (LD), respectively. The angle of local moon normal to the star was small, thus this 
observation sampled the minor axis of the star. A second set of observations obtained 
several months later at KPNO by Ridgway ct al. (1982) with the 4-m telescope in //-band 
(1.67 iim) lead to a diameter of 1.44±0.09 mas (UD) or 1.46 (LD). The angle of approach was 
somewhat steeper in this instance, about 32° (D. Dunham, 2004, private communication). 
It is tempting to suppose that the somewhat higher value in this case is due to the geometry, 
but the observations obviously agree within stated errors. The simple average of all three 
observations (averaging the two from Radick) is 1.42 mas (LD), which agrees well with our 
diameter measurements discussed below. 

A detailed consideration of the response of a long-baseline interferometer to a rotation- 
ally distorted star was first undertaken by Johnston & Wareing (1970) who identified Regulus 

and Altair (a Aql), along with the somewhat fainter star ( Oph, as nearly ideal targets for 
an interferometer with sensitivity beyond that of the Narrabri instrument. They argued that 
the high projected rotational velocities of these stars implied inclinations of their rotational 
axes with respect to the line of sight of nearly 90°. 

Two rapidly rotating stars have now been observed with the new generation of inter- 
ferometers envisioned by Johnston & Wareing (1970). Van Belle et al. (2001) have used the 
Palomar Tcstbed Interferometer (PTl) to measure the oblateness of Altair and to derive a 
Vsini independent of spectroscopy on the basis of a rapidly rotating Roche model to which 
the measured interferometric visibilities were fit. The PTI results led to the determination 
of the star's polar diameter of 3.04±0.07 mas with an equatorial diameter 14% larger. The 
Very Large Telescope Interferometer (VLTI) has been used by Domiciano de Souza et al. 
(2003) to study the oblateness of the Be star Achernar {a Eri). A rotating Roche model 
with gravity darkening was applied to the collective VLTI visibilities following an initial sim- 
plified uniform disk diameter fit to individual visibilities with the conclusion that the Roche 
approximation does not pertain in the case of Achernar. Achernar was found to have a polar 
angular diameter of 1.62±0.06 mas with an extraordinary oblateness ratio of 1.56±0.05. 

The smaller diameter of Regulus in comparison with Altair and Achernar calls for an 
interferometer with baselines longer than those possessed by the PTl and VLTI instruments. 
At the longest baselines of the CHARA Array, access is obtained in the K-hsmd infrared 



-6- 



to visibilities descending to the first null in the regime in which the departures from a 
uniform disk expected for a rapid rotator become increasingly more pronounced. With this 
in mind and inspired by the interest generated from the PTI and VLTI results, Regulus was 
selected as an ideal target for an extensive observing campaign with the CHARA Array. 
The combination of the resulting CHARA interferometric data with a number of constraints 
resulting from spectroscopy of Regulus have permitted us to determine the star's polar and 
equatorial angular diameters and temperatures, inclination angle of the spin axis with respect 
to the line of sight, position angle of the spin axis in the plane of the sky, equatorial rotational 
velocity, and gravity darkening index. 



2. Interferometric Observations 

Regulus was observed on a total of 22 nights between 2004 March 10 and 2004 April 16 
using 10 of the 15 baselines available from the CHARA array. A cahbrator star, HD 83362 
(G8 HI, V — +6.73, K — +4.60), was selected on the basis of its predicted small angular 
diameter and its lack of known stellar companions. This star is located 7?4 from Regulus. No 
suitable cahbrator was found closer in angular separation. A third object, HD 88547 (KO HI, 
V = +5.78, K — +2.97), was chosen as a check star on the basis of having a predicted angular 
diameter comparable to that of Regulus but not expected to show rotational oblatcness due 
to its low V sini = 2.5 + 0.8 km s^^ (Henry et al. 2000). These three objects were observed in 
a sequence so as to provide a series of time-bracketed observations permitting the conversion 
from raw to calibrated visibility for each observation of Regulus and the check star. 

All observations were obtained in the K-hand infrared using a K'-hand filter whose 
effective wavelength was determined to be 2.1501 /im for a star of spectral type B7. This value 
is based upon the atmospheric transmission, the measured filter transmission, the detector 
DQE, and from data provided by the vendors for the various optical surfaces encountered 
along the hght path of each arm of the interferometer. The stellar photon count rate is also 
included, and the increase in effective filter wavelength to 2.1505 //m for the late-type giants 
is considered to be neghgible here. 

Interference fringes were obtained using a classical two-beam interferometer configura- 
tion in which hght is detected emerging from both sides of a beam splitter with a PICNIC 
HgCdTe 256x256 hybrid focal plane array (sensitive to the wavelength range 1.0 to 2.5 //m) 
developed by Rockwell Scientific Company (Sturmann et al. 2002). After open-loop path 
length difference compensation is attained with the optical delay lines, the zero path dif- 
ference position is scanned with a dither mirror in order to produce fringes of a selectable 
frequency, with 150 Hz fringes being utilized for all the data in this analysis. A single dataset 
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consists of some 200 fringe scans preceded and followed by shutter sequences to measure the 
light levels on each input side of the beam splitter. These sequences permit the calculation of 
a visibility adjustment factor to correct for any imbalance between the two detected signals 
(Traub 2000). Such a dataset requires approximately five minutes of observing time. 

2.1. Visibility Measurement 

Methods of extraction of visibihties from these data are discussed in detail in Paper II 

(ten Brummelaar et al. 2004) and generally follow the procedures described by Benson, 
Dyck, & Howell (1995). For the data obtained here, the correlation u, measured from fringe 
amplitude, is treated as the proxy for visibility V. Algorithms were developed for seeking 
the interference fringe in each scan by locating the maximum excursion from mean intensity 
within the scan. Identification of false features in each scan can be a problem at low SNRs, 
and several methods were applied to guard against this. These include the tracking of the 
location of the fringe center, which should wander through successive scans in a continuous 
fashion, and calculating a SNR for each detected fringe based upon the amphtude of the 
peak in the power spectrum for that fringe. When a detected fringe is considered to be 
real, the amplitude of that fringe is determined by mathematical fitting of the fringe, and 
the amplitude of the fit is adopted as the correlation u. This means that the entire fringe 
packet, not just the maximum in the packet, contributes to the measurement of u. The 
sequence of fringe scans in a given dataset is subdivided by time into 15 subsets from which 
15 mean values of u are calculated, and the final value of u and its standard error associated 
with a given dataset are taken as the mean and standard deviation arising from the 15 
subsets. 

The conversion of a correlation measured from fringe amplitude to a visibility is achieved 
by dividing the correlation of the calibrator star obtained by time interpolation between 
cahbrator observations immediately before and after the target star. This is the standard 
practice in optical interferometry to compensate for the time varying transfer function arising 
primarily from atmospheric fluctuations. It assumes a known value for the angular diameter 
of the calibrator whose predicted visibility can then be multiplied by the ratio of the target 
to calibrator correlations to obtain the visibility of the target star. 

Estimating the visibility of the cahbrator star (in this case HD 83362) is a critical 
factor in any interferometric diameter analysis. We have obtained a reliable estimate of 
the diameter of the calibrator by utilizing the relationship existing among the parameters 
angular diameter, effective temperature, and surface flux. Thus, we have calculated an 
estimated angular diameter for the calibrator star HD 83362 by fltting to a spectral energy 
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distribution selected according to spectral type from a template compiled by Pickles (1998) 
to the available broadband photometry, particularly in the near-infrared (Gezari et al. 1993; 
Cutri et al. 2003). The resulting angular diameter for HD 83362 is 0.516±0.032 mas (UD), 
and we adopted this value for the calibration of visibilities for Regulus and the check star 
HD 88547. The error in the estimated diameter for the calibrator star is ±6.2%, and this 
propagates into the measurements for the check star to produce ±0.9% and ±2.5% errors in 
calibrated visibility for baselines of 200 m and 300 m, respectively. For the larger diameter 
of Regulus, these errors are reduced to ±0.5% and ±2.2%. The errors in derived diameters 
resulting from the uncertainty in the diameter of the cahbrator are ±0.8% for the check 
star and ±0.5% for Regulus. The overall errors in these determinations are significantly 
larger and dominated by the random error in visibility measurements. The uncertainty in 
the adopted diameter for the calibrator has a second order, and therefore negligible effect, 
on interferometric determination of relative shape and relative surface brightness variation. 

This approach has led to 69 measurements of calibrated visibilities for Regulus and 40 
such visibilities for HD 88547. Those values are summarized in Tables 1 and 2 and are 
available on the CHARA website in the standard optical interferometry FITS format (Pauls 
et al. 2004). The distributions in the projected baseline plane are shown in Figure 1. These 
data are shown plotted in Figure 2 along with the best fits of uniform disk diameters using 
the standard equation for fitting V{B) to obtain Qud (cf. equation 3.41 of Traub 2000). 
These fits yield Qud,hd88547 — 1.29±0.07 mas and QuD,Reguius — 1.47±0.12 mas. It is 
apparent from the standard deviations of these values and simple inspection of Figure 2 
that, while the diameter fit to the check star is reasonable and in good agreement with the 
diameter predicted from the available photometry, the fit to Regulus is much poorer with 
large systematic residuals. This immediately suggests that the interferometry of Regulus 
cannot be explained in terms of a simple uniform disk model. 



2.2. Modelling Visibilities Geometrically 

The interferometric analyses of the rapid rotators Altair (van Belle et al. 2001) and 

Achernar (Domiciano de Souza et al. 2003) discussed the shapes of these stars on the basis 
of uniform disk diameter values associated with each prior to subjecting their data to 
more sophisticated models. Using this same approach as a starting point, we have selected 
observations made at the longest baselines that arc chistcrcd in several position angle regions 
as indicated in Figure 1. There are 6, 5, and 6 measurements of V respectively in these 
clusters, and they yield UD diameters for Regulus of 1.413±0.024 mas, 1.514±0.028 mas, 
and 1.328±0.028 mas for the three mean position angles 129^1, 181?3, and 251?6, respectively. 
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A similar treatment of observations of the check star HD 88547 results in UD diameters of 
1.233±0.012 mas, 1.260±0.057 mas, and 1.246±0.015 mas at position angles 128?8, 206?8, 
and 253?4. These diameters are shown plotted in Figure 3 which clearly suggests that the 
check star is round while Regulus exhibits an elongated shape whose major axis is oriented 
along a position angle roughly in the range of 150° to 180°. 

We next fitted the visibility data with a series of elliptical models of constant surface 
brightness, which we label as UE fits to distinguish from UD models. These models have the 
advantage that their Fourier transform is just an Airy function, elongated according to the 
axial ratio (Born & Wolf 1999). These fits arc a natural extension to the UD models done 
traditionally for (circular) stellar diameters. Indeed, the fact that UE models arc, like UD 
fits. Airy functions is why one may reasonably take the approach of the previous paragraph 
for a star whose shape is approximated by an ellipsoid. This three-parameter approach yields 
a best-fit UE model with -Rminor/-Rmajor = 0.845 ± 0.029, a mean diameter of 1.42 ± 0.04 mas 
(or -Rminor = 0.651, i?major = 0.771 mas), a position angle of the short axis of the ellipse of 
84?90 ± 2?4 (measured eastward from north), and a reduced chi-square value of 3.41. These 
results indicate that the star is clearly non-spherical in shape, and in the following section 
we develop a physical model for the star so that we can make a realistic comparison between 
the predicted rotationally distorted shape and the observed visibilities. 



3. Spectroscopic Constraints on the Physical Pctrameters of Regulus 

Rapid rotation in stars like Regulus has two immediate consequences for a physical 
description of the star. First, the star will become oblate, and thus we need a geometrical 
relationship for the star's radius as a function of angle from the pole. Secondly, the surface 
temperature distribution will also become a function of co-latitude, usually assumed to be 
cooler at the equator than at the poles. Rehable estimates of the predicted stellar flux, spec- 
tral features, and angular appearance in the sky will then be functions of nine parameters: 
mass M, polar radius Rp, polar effective temperature Tp, equatorial velocity V^, inclination 
angle of the pole to the line of sight i, position angle a of the polar axis from the north 
celestial pole through the east at the epoch of observation, gravity darkening exponent /3 
that defines the amount of equatorial cooling (Collins et al. 1991), distance to the star d, and 
interstellar extinction. In principle one can make a grid search through these parameters to 
find the best-fit of the interferometric visibilities for different baselines projected onto the sky 
(van Belle et al. 2001), but in the case of Regulus there are a number of spectral observations 
available that can help to reduce significantly the probable range in these parameters. Here, 
we discuss how the spectral ffux distribution and the appearance of several key line profiles 



-10- 



can be analyzed to form these additional constraints that we will apply to the analysis of 
the interferometric visibility data. 

We developed a physical model for the star by creating a photospheric surface grid of 
40000 area elements. We assume that the shape of the star is given by the Roche model 
(i.e., point source mass plus rotation) for given equipotential surfaces based on the ratio of 
equatorial to critical rotational velocity V^/V^ (Collins et al. 1991). The effective temperature 
is set at each co- latitude according to the von Zeipel law (von Zeipel 1924), in which the 
temperature varies with the local gravity as T(9) = Tp{g{9)/ Qp)^ and (3 is normally set 
at 0.25 for stars with radiative envelopes (Claret 2003). The model integrates the specific 
intensity from all the visible surface elements according to their projected area, temperature, 
gravity, cosine of the angle of the surface normal to the line of sight (/x), and radial velocity. 
The intensity spectra were calculated over the full range of interest in steps of 2000° in 
temperature, 0.2 dex in logg^, and 0.05 in /i. The spectra were computed using the code 
Synspec (Lanz & Hubeny 2003) and are based on solar abimdance atmospheres models by 
R. L. Kurucz'^ (for a chosen microturbulent velocity of 1 km s~^). The model incorporates 
both limb darkening and gravity darkening and provides reliable line profiles even in cases 
close to critical rotation (Townsend, Owocki, & Howarth 2004). Each model spectrum is 
convolved with an appropriate instrumental broadening function before comparison with 
observed spectra. The code is also used to produce monochromatic images of the star in the 
plane of the sky. 

We begin by considering the full spectral fiux distribution to develop constraints on the 
stellar radius and temperature. We plot in Figure 4 observations of the fiux of Regulus binned 
into wavelength increments of the same size as adopted in the Kurucz model flux spectrum 
shown. The measurements for the extreme ultraviolet are from Morales et al. (2001), and 
the far- and near-ultraviolet spectra are from the archive of the International Ultraviolet 
Explorer (lUE) Satellite (spectra SWP33624 and LWP10929, respectively). The optical 
spectrophotometry is an average of the observations from Alekseeva et al. (1996) and Le 
Borgne et al. (2003). The near-infrared fluxes are from the IR magnitudes given by Bouchet 
et al. (1991) that we transformed to fluxes using the calibration of Cohen et al. (1992). We 
initially flt this distribution using simple flux models for a non-rotating star, and the flt 
illustrated in Figure 4 uses the parameters Teg = 12250 K, \ogg = 3.5, E{B -V) ^ 0.005, 
and a limb-darkened angular diameter of 1.36 ± 0.06 mas. These parameters agree well with 
those from Code et al. (1976) based upon a full flux integration and the angular diameter 
from intensity interferometry. 
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We used the fit of the spectral energy distribution in X-band as the basis of an obser- 
vational constraint on the integrated flux for our rotating star model. The comparison was 
made by calculating the K-hand flux for a given model and then adjusting the flux accord- 
ing to the distance from the Hipparcos parallax measurement of Regulus, tt = 42.09 ± 0.79 
mas or d = 23.8 ± 0.4 pc (Ferryman et al. 1997), and the negligible extinction derived from 
E{B — V) (Fitzpatrick 1999). We require that modeled and observed fluxes agree to better 
than 1% (the error associated with the fit of the absolute observed flux in the if -band). This 
constraint primarily affects the selection of the polar temperature (which sets the flux) and 
polar radius (which together with the adopted distance sets the angular area of the source) . 

The next constraint is provided by the shape of the hydrogen H7 A4340 line profile 
as shown in Figure 5. This observed spectrum is from the spectral library of Valdes et al. 
(2004), and it has resolving power of A/ A A = 4900. This Balmer line grows in equivalent 
width with declining temperature through the B-star range, and it is wider in dwarfs than 
in supergiants because of pressure broadening (linear Stark effect). Thus, fits of the profile 
are set by the adopted model temperature and surface gravity distribution over the visible 
hemisphere of the star. 

Our final constraint is established by the rotational line broadening observed in the 
Mg II A4481 line as shown in Figure 6. This spectrum is the sum of 30 spectra obtained in 
1989 April with the Kitt Peak National Observatory 0.9-m Coudc Feed Telescope, and it has 
a resolving power of A/AA = 12400. This spectral region also includes the weaker metallic 
lines of Ti II A4468, Fe II A4473, and Fe I A4476, which become prominent in cooler A- 
type spectra. All these lines have intrinsically narrow intensity profiles, and the broadening 
we observe is due primarily to Doppler shifts caused by stellar rotation. A fit of the Mg II 
A4481 profile leads directly to the projected equatorial rotational velocity, F sin i, which for a 
given inclination then defines the actual rotational speed and stellar deformation. Note that 
we exclude the nearby He I A4471 line from this analysis, since it is possible that Regulus 
belongs to the class of He- weak stars that account for about one quarter of the B-stars with 
temperatures similar to that of Regulus (Norris 1971). 

We assume throughout our analysis that the optical and if -band fiux originates solely in 
the stellar photosphere. This assumption needs verification because many rapidly rotating B- 
type stars develop large equatorially confined disks (spectroscopically identified as Be stars; 
Porter & Rivinius 2003) and such disks can contribute a large fraction of the if -band flux 
(Stee & Bittar 2001). The appearance of excess IR flux from a disk is always accompanied by 
the development of an emission feature in the Balmer Ha proflle (Stee & Bittar 2001), so we 
have searched for any evidence of Ha emission in spectra obtained contemporaneously with 
the interferometric observations. We show in Figure 7 the average Ha profile in the spectrum 
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of Regulus formed from 11 spectra obtained with the KPNO Coude Feed Telescope on 2004 
October 13 - 16 (less than six months after our most recently collected interferometric data). 
This spectrum has a resolution of i? = A/AA = 10600 and a signal-to-noise ratio of 500 per 
pixel in the continuum. It appears completely free of any excess emission such as that 
shown in Figure 7 for two Be stars observed at the same time as Regulus (HD 22780 and 
HD 210129, both of comparable spectral classification, B7 Vne). We also show in Figure 7 
the photospheric Ha profile of a slowly rotating B-star, HD 179761, which we obtained from 
the spectral atlas of Valdes et al. (2004). This star has a temperature — 13000 K and 
a gravity logg' = 3.5 (Smith & Dworetsky 1993) that are comparable to the surface mean 
values for Regulus. However, it has a very small projected rotational velocity, Vsmi — 15 
km s"^ (Abt, Levato, & Grosso 2002), so we artificially broadened the profile by convolution 
with a rotational broadening function with V sin i = 317 km s^^ (see below) and a linear limb 
darkening coefficient of 0.30 (Wade & Rucinski 1985) in order to match approximately the line 
broadening in the spectrum of Regulus. The good agreement between the spun-up version 
of the Ha profile of HD 179761 and that of Regulus indicates that no excess disk emission 
is present. We have also examined spectra {R — 14000; S/N = 200 - 500) made with the 
CHARA 1.0-m equivalent aperture Multi- Telescope Telescope (MTT) (Barry, Bagnuolo, & 
Riddle 2002) on six nights (1999 January 19, March 17, 2000 February 29, April 16, 20, 29), 
and none of these showed any evidence of emission. Furthermore, another KPNO Coude 
Feed spectrum from 2000 December and 21 HEROS spectra (Stefi, Hummel, & Rivinius 
2000) made between 2001 May and 2002 May (kindly sent to us by Dr. S. Stefi) are free of 
Ha emission. These contemporary observations of no emission are generally consistent with 
the record in the literature. For example, Slettebak & Reynolds (1978) used Regulus as a 
standard emission-free and non-varying star in a survey for Ha variations among Be stars 
during the period 1975 December to 1977 June. There is only one report of marginal Ha 
emission in the spectrum of Regulus in 1981 February (Singh 1982), but the observations at 
hand indicate that there was no disk at the time of the CHARA Array observations, so that 
the K-hand visibilities can be analyzed reliably in terms of a photospheric model alone. 

We constructed a series of stellar models that consistently meet all of these observa- 
tional constraints for a grid of assumed values of stellar inclination i and gravity darkening 
coefficient f3. The method was an iterative approach in which we first established the pro- 
jected rotational broadening from fits of the Mg II A4481 line (dependent only on over 
the inclination range of interest) and then progressively approximated the polar temperature 
and radius to find models that met both the A'-band fiux and H7 profile constraints. Note 
that once the temperature and radius are set, the H7 fit provides us with the gravity (log 51) 
and hence the final parameter, stellar mass. The results of these solutions for the i and P 
grid are listed in Table 3, which includes the fitted projected rotational velocity (column 3), 
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the ratio of equatorial to critical velocity (column 4) , the polar and equatorial radii (columns 
5 and 6), the mass (column 7), and the polar and equatorial temperatures (columns 8 and 
9). 

The first derived quantity is the projected rotational velocity Vsini, which depends 
mainly on the adopted value of gravity exponent /3. With larger (3, the equatorial regions 
are darker and contribute less in the extreme line wings. Consequently, the best fits ac- 
commodate a slightly larger value of Vsini at the nominally accepted value of f3 = 0.25 
compared to the case of no gravity darkening, /? = 0. Fitting errors admit a range of ±3 
km s~^ in the derived V sini value. We find, however, that the fits of Mg II A4481 become 
much worse for models close to critical rotation, which confirms the inclination range of 
60° < i < 90° found by Stoeckley & Buscombe (1987) based on similar profile studies. The 
estimates of stellar radii and mass are almost independent of the choice of i and /3, and it 
is only the range in stellar temperature between the equator and pole that increases with 
/3. Note, however, that the geometric mean of the equatorial and polar temperatures is 
almost the same in all the models, which results from matching the integrated /T-band flux 
(reflecting the average temperature across the visible hemisphere) . The internal errors of the 
fitting procedure (exclusive of any systematic errors associated with the model atmospheres 
and fluxes) lead to errors of ±1% in the radii and ±7% in mass. 



4. Fits of Interferometric Visibility Using Physical Models 

All of the models in Table 3 make acceptable fits of the spectroscopic data, and to 
further discriminate between them we turn to the model predictions about the interferometric 
visibilities. Here we discuss how to predict the visibility patterns associated with the physical 
models and how the observed visibilities can be used to evaluate the remaining parameters, in 
particular the inclination and the gravity darkening exponent. We first explore the goodness 
of fit for the full grid of models shown in Table 3 and we determine what parameter set 
minimizes the reduced chi-square for the whole set and subsets of the visibilities. We then 
compare the results with an independent approach based upon a guided search of the multi- 
parameter space. 

We created X-band surface intensity images of the stellar disk for each case in Table 3 
and made a two-dimensional Fourier transform of the image for comparison with the inter- 
ferometric visibility (Domiciano de Souza et al. 2002). The modulus of the Fourier transform 
relative to that at frequency zero is then directly comparable to the observed visibility at 
a particular coordinate in the {u, v) spatial frequency plane (defined by the ratio of the 
projected baseline to the effective wavelength of the filter). 



-14- 



The predicted variation in visibility as a function of baseline is shown in Figure 8 for three 
selected models and for two orientations: the top set of three plotted lines corresponds to 
an interferometric baseline parallel to the polar (minor) axis while the lower set corresponds 
to the equatorial (major) axis. The CHARA Array observations correspond to baselines 
between 188 and 328 m, so we are sampling the main lobe of visibility. Model A corresponds 
to the case with i = 90° and (3 = 0.25 (see Table 3). The visibility drops much more 
rapidly in the direction along the major axis compared to the minor axis, and consequently 
the variations in interferometric visibility in different directions in the {u, v) plane set the 
orientation of the star in the sky. Model B is also based on the nominal gravity darkening 
exponent value oi P — 0.25, but the inclination is i = 70° in this case. There are subtle 
differences in the pole to equator axis visibilities from model A that result from the greater 
rotational distortion and smaller polar radius in model B. Finally, model C shows the i = 90° 
case but with no gravity darkening {j3 = 0). Now the equatorial flux is easily seen out to 
the extreme approaching and receding limbs, and the star appears larger in this dimension 
(resulting in lower visibility). These models show that although the differences due to stellar 
orientation in the sky are seen at all baselines, the smaller differences related to stellar 
inclination and gravity darkening are most apparent at longer baselines. Indeed, future 
observations with even longer baselines (or shorter wavelengths) that reach into the second 
lobe would strongly discriminate between these models. 

We made a two-parameter fit to match the predicted and observed visibilities of Regulus. 
The first parameter is the position angle of the minor (rotation) axis in the sky a. The second 
parameter is an angular scahng factor, the ratio of best fit distance to the adopted distance 
from the Hipparcos parallax. In principle, the combination of the i^'-band flux and adopted 
temperature model leads to a predicted angular size (or stellar radius for a given distance), 
but the observed errors in absolute flux and distance will lead to a predicted angular size that 
may be marginally different from the actual angular size. The use of the scaling parameter 
provides a simple way to both check the prediction and adjust the spatial frequency scaling to 
best fit the observed visibihties. The results of these fits are given in Table 3 which lists the 
minor axis angular radius (column 10), major axis angular radius (column 11), the position 
angle a (column 12), the scaling factor d/d{Hipparcos) (column 13), and the reduced chi- 
square of the flt xt (column 14; for N — 69 measurements and z/ = 2 fitting parameters). We 
estimate that the errors in position angle are ±2? 8 and in scaling factor are ±1% (based upon 
the increase in xl with changes in these parameters). We see that the scaling factors are 
generally less than 2% different from unity, well within the error budget from the comparable 
errors in the estimates of the i^-band flux, parallax, and effective wavelength of the CHARA 
K' filter. 

The best overall fit is found for the rotational model with i — 90° and j3 — 0.25, and 
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we show the predicted image of the stellar disk in the sky and its associated visibility in the 
(u, v) plane for this model in Figure 9. The upper part of the right panel shows the visibility 
as a gray-scale intensity and the polar dotted line (along the longer dimension in the 

visibility distribution). Black squares mark the positions in the {u,v) plane corresponding 
to the CHARA interferometric observations. The lower panel shows the point symmetric 
{u, v) plane, and there each observation is assigned a gray-scale intensity based upon the 
normalized residual, (Kbs ~ KaO/c^y- Note that the poor fits are the most apparent in this 
representation, but that there is no evidence of any systematic trends in the residuals with 
position in the {u, v) plane. The relative close proximity in the {u, v) plane of measurements 
that are above and below the predicted visibility suggests that the outstanding differences 
are due to measurement error rather than deficiencies in the models. The fact that the 
minimum chi-square is xl = 3.35 rather than the expected value of unity indicates that our 
internal estimates of error in visibility (from the scatter within a sequence of fringe scans) 
may not adequately represent the total error budget. However, it is well known that the 
temporal power spectrum of atmospheric fiuctuations has a steep slope. Since the calibration 
cycle has a much lower frequency than the scan to scan cycle, it is natural that the noise 
estimate from the latter underestimates the former. 

We showed above in Figure 8 that the interferometric visibilities at all baselines provide 
a strong constraint on the position angle a, and we illustrate this striking dependence in 
Figure 10. This series of panels shows the normalized residuals from the model predictions 
for the i — 90° and P — 0.25 case as a function of basehne. The sky orientation of the 
K-hand image of the star is shown to the right of each panel. The residuals for all baselines 
clearly obtain a well defined minimum for a = 85?5 {third panel from top). Figure 11 shows 
the reduced chi-square as a function of a for both the entire set of measurements {solid line) 
and the long baseline measurements only (-Bproj > 270 m; dotted line). Both sets yield a 
consistent minimum, and the long baseline data are particularly sensitive to the position 
angle orientation. 

The visibility constraints on the inclination and gravity darkening exponent are less 
pronounced but still of great interest. We show in Figure 12 the reduced chi-square as a 
function of the gravity darkening exponent /? for a series of i = 90° model fits. The best fit 
occurs at P — 0.25 for the full set of observations, and this is also the value derived from 
gravity darkening studies of B-stars in eclipsing binary stars (Claret 2003). The formal la 
error hmit yields an acceptable range from (3 — 0.12 to 0.34, but P — (no gravity darkening) 
can only be included if we extend the range to the 99% confidence level. However, recall from 
Figure 8 that most of the sensitivity to gravity darkening is only found at longer baselines 
and especially at those along the polar axis. Thus, we also show in Figure 12 the value of 
reduced chi-square for the i = 90° solutions in two subsets: measurements with baselines 
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greater than 270 m (31 points) and those with a position angle close to the polar axis 
orientation (6 points from MJD 53103.2 - 53104.2; see Table 1). Both sets (especially the 
latter, polar one) show much larger excursions and indicate (3 values close to 0.25. Note that 
we did not try models with j3 > 0.35 since such models have equatorial temperatures that are 
cooler than the lower limit of our flux grid. The most sensitive polar subset helps establish 
more practical error limits, and we adopt f3 = 0.25 ± 0.11 as our best overall estimate. This 
analysis demonstrates that the inclusion of gravity darkening is key to fitting the CHARA 
data in the most sensitive part of the visibility curve. 

The greater sensitivity of the long baseline and polar data is also apparent in the flts 
with differing inclination angle. We show in Figure 13 the run of reduced chi-squarc with 
inclination for models with /3 = 0.25. A shallow minimum at i — 90° is found using all the 
data, but the same minimum is much more convincingly indicated in fits of the long baseline 
and polar data alone. The la error hmit of the latter sample admits a range of 78° — 90°. 
Hutchings & Stoeckley (1977) also derived an inclination of i = 90° for Regulus through a 
comparison of ultraviolet and optical photospheric line widths. 

It is interesting to compare the rotation axis vector direction, determined here by com- 
bination of the angles i and a, with the space velocity vector direction. Hipparcos proper 
motion and parallax results (Ferryman ct al. 1997) combined with a new radial velocity 
value of V = +7.4 ± 2.0 km s~^, derived from the analysis described above of the Mg II 
A4481 line, provide the necessary input for making this comparison. The angular values 
representing the space velocity that correspond to our measurements of i = 90° ± 15° and 
a = 85?5 ± 2?8 are ig = 75?2 ± 4?0 and as = 271? 1 ± 0?1. Because a is inherently ambiguous 
by 180°, a-as^ -5?6 ± 2?8 and i - 14?8 ± 15?5. This suggests that Regulus is moving 
very nearly pole-on through space. 

We have checked these models through a comparison with an independent code based 
upon the scheme first outlined by van Belle et al. (2001), using a new generation of that code 
that conducts its comparisons in Fourier space rather than in image space. This method 
uses a different numerical realization of the geometry of the star, but one that is based on 
the same Roche approximation adopted above. There are six key parameters involved in 
this method (polar radius, i, a, j3, rotation speed, and, albeit with low sensitivity, polar 
temperature) that define the projection of the stellar surface on the sky. This sky projection 
was then run through a two-dimensional Fourier transform, and a reduced chi-square was 
calculated from a comparison of the observed and model values for V"^. A multi-dimensional 
optimization code was then utilized to derive the best solution, a process that took typically 
500 iterations (Press et al. 1992). The results are compared with those from fitting a uniform 
ellipsoid (§2.2) and from the spectroscopically constrained fits (Table 3) in Table 4. We 



-17- 



see that all three methods agree on the orientation of the polar axis in the sky, and both 
the grid search and spectroscopically constrained approaches find the same results for the 
angular sizes and inclination. The only discrepancy concerns the value of j3, and the lower 
value derived from the grid search scheme is probably the result of the insensitivity to this 
parameter of the majority of the measurements in the whole sample. We place more reliance 
on the results of the sensitive polar sample in the spectroscopically constrained fit (Fig. 12). 

In summary, we find that infrared interferometric measurements of Regulus over a wide 
range of position angle are consistent with parameters derived from spectroscopic criteria, 
and determine additional parameters which arc not available from spectroscopy, particularly 
including the position angle of the rotation axis on the sky a, the inclination of that axis 
to the line of sight i, and the gravity darkening coefficient (5. Our adopted results and their 
errors are summarized in Table 5. Our physical models for this rapidly rotating star indicate 
that Regulus has an equatorial radius that is 32% larger than the polar radius. Its rotation 
period is 15.9 hours, which corresponds to an equatorial rotation speed that is 86% of the 
critical break-up velocity. Fits of the CHARA observations require the presence of significant 
gravity darkening with an equatorial temperature that is only 67% of the polar temperature. 



5. Conclusions 

We have shown in a series of increasingly complex models, ultimately tied to physical 
parameters strongly constrained by spectroscopy, that Regulus exhibits features expected 
for a rapidly rotating star of its spectral type, namely oblateness and gravity darkening. 
Geometric fits to fringe visibility, first with discrete position angle determinations of uni- 
form disk diameters and then with an ellipsoidal model, clearly show the marked rotational 
oblateness of the star. When we couple the visibilities to models that incorporate parameters 
to which high resolution spectroscopy is sensitive, we find mutual consistency between the 
interferometric and spectroscopic results. Indeed, we believe that the combination of these 
complementary astrophysical probes - interferometry and spectroscopy - provides the best 
means for exploiting new high spatial resolution measurements from such instruments as the 
CHARA Array. 

Our observations provide the first interferometric evidence of gravity darkening in 
rapidly rotating stars. Furthermore, the CHARA Array results offer the first claim of grav- 
ity darkening in a star that is not a known member of an eclipsing binary system (Claret 
2003). The agreement between the measurement of the angular diameter from the CHARA 
visibilities and that based upon the i^T-band fiux provides an independent verification of the 
infrared flux method for estimating angular diameters flrst introduced by Blackwell & Shallis 
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(1977), and it indicates that the B-star fluxes predicted by hne blanketed, LTE atmospheres 
models by R. L. Kurucz agree with the observed angular diameter and infrared flux within 
the observed errors. 

New interferometric observations of this kind offer the means to determine how close the 
most rapidly rotating stars are to their critical rotation speeds, and this may help solve the 
longstanding problem of the nature of mass loss in the rapidly rotating Be stars, for example 
(Porter & Rivinius 2003). These observations may also finally allow us to test models of 
interior structure and evolution for massive rotating stars (Endal & Sofia 1979; Hcger & 
Langcr 2000; Mcynct & Macdcr 2000). We show in Figure 14 the position of Rcgulus in 
the theoretical Hertzsprung-Russcll diagram based upon the surface integration of crT^ to 
estimate the luminosity and upon an average temperature derived from the luminosity and 
surface area of our best-flt model. We also show the zero-age to terminal-age main sequence 
evolutionary paths for non-rotating 3 and AMq stars based upon the work of Schaller et 
al. (1992). The shaded region represents the main sequence evolutionary parameters for a 
star with the mass and errors in mass we derive from the joint spectroscopic-interferometric 
analysis. The position of Regulus in this diagram shows that the star is overluminous for 
its mass as predicted for rotating stars that evolve to higher luminosity as fresh H is mixed 
into their convective cores (Heger & Langer 2000; Meynet & Maeder 2000). If Regulus is 
placed on an isochrone without consideration of this excess luminosity, the age of the star 
is overestimated. This is the likely explanation for the apparent discordance in age that has 
been noted for the a Leo A and B components (Gerbaldi et al. 2001). The new era of long 
baseline interferometry offers us the means to probe the evolution of rotating stars as has 
never before been possible. 
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Fig. 1. — The baseline coverages for Regulus {open diamonds) and the check star HD 88547 
{plus signs) are shown in the {u, v) plane above. The lighter shaded regions in the lower 
half indicate those measurements of Regulus used for a position angle dependent estimate 
of diameter (see Fig. 3), while the darker shaded regions in the upper half show the same 
for measurements of the check star. 
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Fig. 2. — The best-fit uniform disk diameter fits to tlie calibrated visibilities are shown for 
the check star HD 88547 {top panel) and Regulus {bottom panel). The data for Regulus are 
clearly poorly fit by a uniform disk model. 
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Fig. 3. — Uniform disk diameters calculated from the longest baseline data in the vicinity 
of three position angles indicated in Fig. 1 are shown for the check star HD 88547 {lower 
values) and Regulus {upper values). The dotted lines show the respective mean values of 
these 3 sets in each case. These results imply that the check star is round while Regulus is 
not. 
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WAVELENGTH (microns) 

Fig. 4. — The full spectral flux distribution of Regulus. Smaller plus signs indicate the short 
wavelength observations which are binned to the same resolution as the model distribution 
{solid line) while the larger plus signs indicate fluxes from infrared magnitudes. 
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Fig. 5. — The observed {plus signs) and model {solid line) profiles for H7 in the spectrum of 
Regulus. The plot corresponds to a model with a gravity darkening exponent P — 0.25 and 
an inclination i — 90°. The lower plot shows the residuals from the fit (note that the He I 
A4387 line that was not included in the calculation) . 
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Fig. 6. — The observed {plus signs) and model {solid line) profiles for Mg II A4481 in the 
spectrum of Regulus. The plot corresponds to a model with i = 90° and (3 = 0.25. Residuals 
from the fit are shown below. Note that the He I A4471 line and other weaker lines were not 
included in the fit of projected rotational velocity (which was based on the interval between 
the two vertical dashed lines.) 
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Fig. 7. — A plot of the Ha profile in the spectrum of Regulus {bottom; HD 87901). This 
is an average of eleven spectra made with the KPNO Coude Feed Telescope from 2004 
October 13 - 16. There is no evidence of the kind of disk Ha emission that is observed in 
rapidly rotating Be stars such as HD 22780 {second from top; B7 Vne) and HD 210129 {top; 
B7 Vne). Also shown is the photospheric Ha profile of the slowly rotating star HD 179761 
{third from top; B8 II-IH), which we artificially broadened to match the rotational broadening 
of Regulus. The close agreement between the photospheric line of HD 179761 and that of 
Regulus indicates that no disk emission is present. 
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Fig. 8. — Predicted visibility variations with baseline for three rotation models and two sky 
orientations. The upper group corresponds to a baseline parallel to the minor (rotational) 
axis in the sky while the lower group corresponds to a baseline parallel to the major axis. The 
visibility curves are shown for the cases oil = 90°, i3 = 0.25 (A; solid lines), i = 70°, f3 = 0.25 
(B; dotted lines), and i = 90°, f3 = (C; dashed lines). The spatial frequency for a given 
baseline is shown for a filter effective wavelength of 2.1501 /xm. The shaded region indicates 
the basehne range of the CHARA Array observations. 
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Fig. 9. — A A'-band image of the star in the sky (left) and its associated Fourier transform 
visibihty pattern in the {u,v) plane (right). In both cases north is at the top and east is to 
the left. The dotted, black line indicates the direction of the rotational axis for this i — 90°, 
P — 0.25, and a — 85?5 model. The upper panel of the visibihty figure (right) shows a 
grayscale representation of the visibihty and the positions of the CHARA measurements 
(black squares). The lower panel shows the normalized residuals from the fit as a gray scale 
intensity square against a gray background in a point symmetric representation of the (u, v) 
plane. The legend at lower left shows the intensities corresponding to normalized residuals 
from —5 (black) to +5 (white). Note that the best fit points appear gray and merge with 
the background. 
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Fig. 10. — Normalized visibility residuals as a function of baseline. Each panel shows the 
residuals for the model star with i = 90°, f3 = 0.25, and a position angle a as indicated (and 
illustrated at right). The residuals are clearly minimized at the best fit value of a = 85?5 
{third panel from top). Plus signs indicate measurements in the {u,v) plane within 30° of 
the rotation axis (6 points), diamonds indicate those within 30° of the equator (40 points), 
and asterisks indicate the others at intermediate angles (23 points). 
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Fig. 11. — A plot of the reduced chi-square xl of the visibihty fits as a function of position 
angle a (for i = 90° and f3 = 0.25). The solid line shows the reduced chi-square for whole 
sample while the dotted line shows the same for the long baseline data only. 
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Fig. 12. — A plot of the reduced chi-square xl of the visibihty fits as a function of gravity 
darkening exponent f3 (for i = 90°). The sohd hne shows the reduced chi-square for whole 
sample (69 points). The dotted line shows the reduced chi-square of the same fits for the 
long baseline data only (31 points) while the dashed line shows the same for those long 
baseline data with a position angle near the orientation of the polar axis (the 6 points most 
sensitive to the selection of /3). All these samples indicate a gravity darkening exponent near 
the predicted value oi P — 0.25. 
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Fig. 13. — A plot of the reduced chi-square xl of the visibihty fits as a function of the 
rotation axis inchnation angle i (for (3 = 0.25). The different lines correspond to the same 
samples as shown in Fig. 12. 
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Fig. 14. — The Hertzsprung- Russell diagram for Regulus. The two sohd hnes show the track 
from zero age to terminal age main sequence {left to right) for non-rotating stars with initial 
masses of 3 and AMq (Schaller et al. 1992), while the dotted line and surrounding shaded 
area show the predicted region for a non-rotating star with the derived mass of Regulus and 
its associated error. The single point at L/Lq = 347 and < T^s >= 12901 K (from a surface 
integration of our numerical model) is well above the predicted position, indicating that the 
star is overluminous for its mass in comparison to models for non-rotating stars. 
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Table 1. Interferometric Measurements of Regulus (HD 87901) 



Projected Baseline 

Baseline Position Angle Calibrated 
MJD (m) (deg) Visibility 



53074, 


.188 


188, 


.76 


181, 


,.33 


0. 


,520 




0.049 


53074, 


.200 


188, 


.77 


178, 


,49 


0. 


,508 


± 


0.069 


53074, 


.215 


189, 


.31 


174. 


.62 


0. 


,535 


± 


0.073 


53074 


.227 


190 


.11 


171. 


.71 


0. 


,613 


± 


0.041 


53074, 


.238 


191, 


.17 


168, 


.98 


0. 


,576 


± 


0.052 


53074, 


.250 


192, 


.51 


166, 


,29 


0. 


,556 


± 


0.071 


53074, 


.263 


194, 


.26 


163, 


,41 


0. 


,476 


± 


0.084 


53078, 


.237 


222, 


.54 


140. 


.50 


0. 


,486 


± 


0.032 


53078, 


.262 


232, 


.97 


135. 


.70 


0. 


.539 


± 


0.039 


53078, 


.280 


239, 


.36 


132. 


.82 


0. 


,528 


± 


0.041 


53078 


.297 


244, 


.18 


130, 


.51 


0. 


,483 


± 


0.033 


53079, 


.207 


210, 


.30 


146, 


,95 


0, 


,546 


± 


0.063 


53079, 


.222 


216, 


.82 


143, 


,36 


0. 


,461 




0.069 


53079, 


.235 


222, 


.84 


140, 


,36 


0. 


,427 


± 


0.059 


53079 


.269 


236, 


.42 


134. 


.15 


0. 


,490 


± 


0.044 


53080, 


.250 


233, 


.17 


202. 


.38 


0. 


.405 


± 


0.045 


53080, 


.309 


223, 


.41 


189. 


.25 


0. 


,334 


± 


0.040 


53080 


.351 


221, 


.40 


178, 


.92 


0. 


,400 


± 


0.038 


53080, 


.366 


221, 


.95 


175, 


,06 


0. 


,433 


± 


0.036 


53080, 


.382 


223, 


.23 


171, 


,17 


0. 


,440 




0.039 


53081, 


.167 


246, 


.64 


214, 


,46 


0. 


,289 


± 


0.031 


53081, 


.181 


245, 


.06 


212, 


.86 


0. 


,310 


± 


0.030 


53081 


.196 


242, 


.83 


210. 


.87 


0. 


,365 


± 


0.036 


53081 


.228 


237, 


.11 


206, 


.00 


0. 


,563 


± 


0.054 


53081, 


.243 


234, 


.04 


203, 


.22 


0. 


,553 


± 


0.048 


53081 


.2.59 


230 


.92 


200, 


,10 


0. 


,524 


± 


0.0.54 


53081, 


.288 


225, 


.89 


193, 


.79 


0. 


,560 


± 


0.087 


53081, 


.307 


223, 


.35 


189. 


.13 


0. 


,509 


± 


0.064 


53081, 


.339 


221, 


.40 


181. 


.09 


0. 


,376 


± 


0.029 


53081 


.354 


221, 


.53 


177, 


.40 


0. 


,387 


± 


0.025 


53081, 


.368 


222, 


.29 


173, 


,79 


0, 


,420 


± 


0.031 


53088, 


.224 


313, 


.63 


207, 


,81 


0. 


148 


± 


0.012 


53088, 


.243 


307, 


.88 


204, 


,22 


0. 


172 


± 


0.016 


53088, 


.262 


302, 


.17 


200. 


.16 


0. 


,178 


± 


0.012 


53088, 


.281 


296, 


.95 


195. 


.59 


0. 


.200 


± 


0.011 


53088, 


.299 


293, 


.18 


191. 


.18 


0. 


,203 


± 


0.013 


53088 


.314 


290 


.95 


187, 


.43 


0. 


,209 


± 


0.012 


53088 


.329 


289, 


.54 


183, 


,33 


0. 


,205 


± 


0.013 


53088, 


.340 


289, 


.19 


180, 


,37 


0. 


,208 


± 


0.013 


53088, 


.351 


289, 


.39 


177, 


,45 


0. 


198 


± 


0.017 


53092, 


.231 


307, 


.94 


204. 


.26 


0. 


,183 


± 


0.021 


53092, 


.250 


302, 


.49 


200. 


.41 


0. 


.191 


± 


0.016 


53093, 


.148 


328, 


.43 


216. 


.70 


0. 


.115 


± 


0.012 


53093, 


.163 


325 


.90 


214. 


.97 


0. 


.112 


± 


0.013 
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Table 1 — Continued 



Projected Baseline 
Baseline Position Angle Calibrated 
MJD (m) (deg) Visibility 



53093.179 


322.24 


212.79 


0.110 


± 


0.007 


53093.199 


316.93 


209.74 


0.131 


± 


0.010 


53093.223 


309.63 


205.36 


0.146 




0.014 


53093.239 


304.84 


202.15 


0.160 


± 


0.011 


53093.256 


300.01 


198.40 


0.151 


± 


0.018 


53093.335 


289.30 


178.12 


0.137 


± 


0.012 


53095.249 


243.70 


130.75 


0.433 


± 


0.041 


53103.200 


312.56 


254.45 


0.263 


± 


0.026 


53103.219 


307.95 


252.80 


0.234 


± 


0.020 


53103.239 


298.97 


250.74 


0.310 


± 


0.021 


53103.259 


285.76 


248.24 


0.318 


± 


0.044 


53104.216 


307.92 


252.79 


0.262 


± 


0.031 


53104.236 


299.02 


250.75 


0.281 


± 


0.033 


53105.261 


278.11 


130.25 


0.340 


± 


0.044 


53105.279 


278.34 


128.86 


0.320 




0.031 


53106.282 


277.96 


128.48 


0.276 


± 


0.028 


53107.290 


276.77 


127.91 


0.259 


± 


0.026 


53108.255 


278.23 


130.09 


0.284 


± 


0.038 


53108.274 


278.14 


128.63 


0.310 


± 


0.045 


53111.265 


260.88 


187.19 


0.227 


± 


0.022 


53111.281 


259.59 


182.72 


0.189 


± 


0.018 


53111.291 


259.37 


180.40 


0.255 


± 


0.017 


53111.300 


259.61 


177.17 


0.213 


± 


0.021 


53111.314 


260.72 


173.20 


0.212 


± 


0.016 


53111.324 


262.02 


170.49 


0.207 


± 


0.017 



-41 - 



Table 2. Interferometric Measurements of the Check Star (HD 88547) 



Projected Baseline 
Baseline Orientation Calibrated 



MJD 


(m) 


[ae 


g) 




Visibility 


53074, 


.369 


205, 


,43 


145, 


.59 


0. 


.560 


± 


0.068 


53074. 


.384 


207, 


,80 


144, 


.20 


0. 


.519 


± 


0.074 


53078. 


.246 


214, 


,24 


136, 


.51 





.563 


± 


0.033 


53078. 


.267 


225, 


,06 


133, 


.10 





.551 


± 


0.039 


53078. 


.286 


233, 


,75 


130, 


.61 


0, 


.532 




0.060 


53079. 


,222 


201, 


,50 


141, 


.20 


0, 


.666 




0.097 


53079. 


,250 


217, 


,65 


135, 


,41 


0, 


.519 


± 


0.068 


53080. 


,267 


217, 


,75 


200, 


,73 


0, 


.522 


± 


0.050 


53080. 


,331 


206, 


,12 


185, 


,08 





.585 


± 


0.049 


53081. 


,181 


240, 


,14 


214, 


,04 


0, 


.512 


± 


0.048 


53081. 


,208 


233 


,51 


210, 


,73 





.455 


± 


0.045 


53081. 


,243 


223, 


,58 


204, 


,98 


0, 


.591 


± 


0.070 


53081. 


,287 


212, 


,36 


195, 


,60 


0, 


.548 


± 


0.081 


53081. 


,322 


206, 


,53 


186, 


,73 


0, 


.574 




0.064 


53081. 


,354 


205, 


,36 


178, 


,02 





.528 


± 


0.033 


53088. 


,230 


298, 


,34 


208, 


,59 





.300 


± 


0.020 


53088. 


,267 


282, 


,95 


200, 


,93 


0, 


.309 




0.025 


53088. 


,287 


275, 


,86 


195, 


,99 


0, 


.364 




0.019 


53093. 


,153 


322, 


,30 


217, 


,29 


0, 


.279 




0.028 


53093. 


,168 


317, 


,64 


215, 


,72 


0, 


.318 


± 


0.027 


53093. 


,187 


310, 


,57 


213, 


,29 


0, 


.319 


± 


0.026 


53093. 


,211 


300, 


,64 


209, 


,54 





.324 


± 


0.022 


53093. 


,228 


293, 


,42 


206, 


,43 





.401 


± 


0.024 


53093. 


,246 


286, 


,07 


202, 


,73 


0, 


.439 


± 


0.034 


53093. 


,261 


280, 


,27 


199, 


,24 


0, 


.404 


± 


0.039 


53095. 


,254 


239, 


,28 


129, 


,06 


0, 


.474 




0.060 


53103. 


,207 


310, 


,56 


254, 


,96 


0, 


.292 


± 


0.031 


53103. 


,224 


304, 


,43 


254, 


,14 


0, 


.332 


± 


0.036 


53103. 


,245 


292, 


,83 


252, 


,89 


0, 


.374 


± 


0.025 


53103. 


,264 


278, 


,05 


251, 


,37 





.409 


± 


0.050 


53104. 


,222 


304, 


,38 


254, 


,14 


0, 


.351 




0.028 


53104. 


,241 


293, 


,38 


252. 


,94 


0, 


.369 




0.031 


53105. 


,267 


276, 


,40 


129. 


,57 


0, 


.413 


± 


0.051 


53105. 


,285 


278, 


,32 


128. 


,73 


0, 


.440 


± 


0.043 


53106. 


,287 


278, 


,48 


128. 


,55 





.424 


± 


0.050 


53107. 


,296 


278, 


,19 


128. 


,27 


0, 


.427 


± 


0.050 


53108. 


,262 


277, 


,03 


129. 


,35 


0, 


.421 


± 


0.042 


53108. 


,283 


278, 


,49 


128. 


,52 


0, 


.421 




0.056 


53111. 


,271 


238, 


,95 


186. 


,96 


0, 


.522 


± 


0.042 


53111. 


,330 


241, 


,57 


168. 


,73 


0, 


.518 


± 


0.041 



Table 3. Parameters of Stellar Models for Regulus 



From Spectroscopy From Interferometry 

i V sin i Rp Re ^minor ^major ^ 

(deg) P (kms-i) Ve/Vc {Re) (Rq) (Mq) (K) (K) (mas) (mas) (deg) d/d{Hipparcos) xl 



90 


0.00 


305 


0.82 


3.15 


4.05 


3.45 


12280 


12280 


0.624 


0.802 


87.4 


0.988 


3.95 


90 


0.05 


307 


0.82 


3.15 


4.06 


3.45 


12825 


11968 


0.624 


0.806 


87.2 


0.987 


3.76 


90 


0.10 


309 


0.83 


3.14 


4.07 


3.45 


13400 


11643 


0.626 


0.810 


86.8 


0.984 


3.58 


90 


0.15 


310 


0.83 


3.14 


4.08 


3.46 


14006 


11309 


0.625 


0.811 


86.1 


0.984 


3.45 


90 


0.20 


314 


0.84 


3.14 


4.12 


3.43 


14660 


10852 


0.624 


0.816 


86.0 


0.987 


3.38 


90 


0.25 


317 


0.86 


3.14 


4.16 


3.39 


15400 


10314 


0.623 


0.825 


85.5 


0.988 


3.35 


90 


0.30 


323 


0.87 


3.15 


4.22 


3.39 


16235 


9677 


0.619 


0.830 


85.3 


0.994 


3.39 


90 


0.35 


326 


0.88 


3.14 


4.24 


3.39 


17120 


9155 


0.620 


0.836 


84.6 


0.993 


3.55 


80 


0.00 


305 


0.83 


3.13 


4.05 


3.45 


12280 


12280 


0.626 


0.804 


87.6 


0.987 


3.97 


80 


0.05 


307 


0.83 


3.12 


4.06 


3.45 


12836 


11943 


0.627 


0.807 


87.2 


0.985 


3.78 


80 


0.10 


309 


0.84 


3.12 


4.07 


3.45 


13420 


11589 


0.626 


0.809 


86.8 


0.985 


3.61 


80 


0.15 


310 


0.84 


3.12 


4.08 


3.45 


14035 


11220 


0.628 


0.814 


86.6 


0.982 


3.48 


80 


0.20 


314 


0.85 


3.12 


4.12 


3.44 


14686 


10740 


0.626 


0.818 


86.1 


0.985 


3.40 


80 


0.25 


317 


0.87 


3.12 


4.16 


3.41 


15435 


10191 


0.625 


0.825 


85.9 


0.987 


3.36 


80 


0.30 


323 


0.88 


3.12 


4.22 


3.40 


16280 


9486 


0.625 


0.835 


85.7 


0.988 


3.36 


75 


0.00 


305 


0.84 


3.10 


4.05 


3.46 


12280 


12280 


0.626 


0.803 


87.7 


0.988 


3.98 


75 


0.05 


307 


0.85 


3.10 


4.07 


3.46 


12845 


11905 


0.629 


0.808 


87.3 


0.985 


3.80 


75 


0.10 


309 


0.85 


3.10 


4.08 


3.46 


13442 


11516 


0.629 


0.810 


86.9 


0.985 


3.64 


75 


0.15 


310 


0.85 


3.10 


4.09 


3.46 


14068 


11109 


0.630 


0.814 


86.7 


0.983 


3.51 


75 


0.20 


314 


0.86 


3.10 


4.12 


3.44 


14740 


10581 


0.630 


0.821 


86.5 


0.983 


3.43 


75 


0.25 


317 


0.88 


3.10 
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Table 4. Parameter Estimates from Different Models 



Uniform Grid Spectroscopically 

Parameter Ellipsoid Search Constrained 



No. parameters 3 6 4 

i^minor (mas) . . 0.651 ± 0.016 0.627 ± 0.015 0.623 ± 0.006 

i?major (mas) . . 0.771 ± 0.032 0.810 ± 0.023 0.825 ± 0.008 

a (deg) 84.9 ±2.4 83.7 ±2.5 85.5 ± 2.8 

i (deg) • ■ • 851^ 90t% 

(3 ••• 0.13 ±0.05 0.25 ±0.11 

Xl 3.41 3.73 3.35 
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Table 5. Final Parameters for Regulus 



Parameter 



Value 



©minor (mas) .. . 1.25 ± 0.02 
©major (mas) ... 1.65 ± 0.02 

a (deg) 85.5 ± 2.8 

^(deg) 90;?5 

(3 0.25 ±0.11 

V'sini (km s-i) 317 ±3 

VjVc 0.86 ±0.03 

Rp (Rq) 3.14 ±0.06 

Re (rJ) 4.16 ± 0.08 

M (Mq) 3.4 ± 0.2 

Tp (K) 15400 ± 1400 

Te (K) 10314 ± 1000 

< T > (K) .... 12901 ± 500 

L (Lq) 347 ±36 

d (pc) 23.5 ± 0.4 

Av (mag) 0.016 ± 0.010 



